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Tau polarization effects in the CNGS v T appearance experiments* 
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We studied r polarization effects on the decay distributions of r produced in the CNGS v T appearance ex- 
periments. We show that energy and angular distributions for the decay products in the laboratory frame are 
significantly affected by the r polarization. Rather strong azimuthal asymmetry about the r momentum axis is 
predicted, which may have observable consequences in experiments even with small statistics. 



1. INTRODUCTION 

An appearance search for — » v T oscillation is 
expected to be performed in the CNGS (CERN 
Neutrino to Gran Sasso) project [I], which will 
start taking data from the year 2006. The project 
with ICARUS and OPERA detectors aims to de- 
tect r-leptons produced via v T J s which originate 
from VaS through oscillation, propagating from 
CERN to Gran Sasso in Italy at a distance of 
732 km. ICARUS observes r decay kinematically 
with liquid argon time projection chamber [2], 
while OPERA does topologically with emulsion 
cloud chamber [3J. 

Because the r is identified by observing its de- 
cay products, and because the decay distribution 
depends crucially on its polarization, it is impor- 
tant to consider the spin polarization of r in ad- 
dition to its production cross section. 

In the last report in the proceedings of 
Nulnt04 [4] , we discussed spin polarization of t~ 
produced in neutrino-nucleon scattering at a fixed 
neutrino energy. Here, we extend to decay dis- 
tributions from produced r~, especially for the 
CNGS experiments. We consider r~ production 
in the neutrino-nucleon scattering and its subse- 
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qucnt decays, for the decay modes r — > 7r v t 
and t~ — ► t~i>(y r (£ — e, fi). 

2. POLARIZATION OF PRODUCED t 

Let us start with a brief summary of spin polar- 
ization of r~ produced via neutrino: (i) at a fixed 
neutrino energy [5] ; (ii) in the CNGS experiments 
with their neutrino flux [B]. 

We consider t~ production by charged current 
(CC) reactions off a nucleon target: 



u T (k)+N(p)^T-(k')+X(p'). 



(1) 



For the hadronic final states X, we consider 
three subprocesses; the quasi-elastic scattering 
(QE), the A resonance production (RES) and the 
deep inelastic scattering (DIS) processes. The 
four-momenta are parametrized in the laboratory 
frame as 



k = (E v , 0, 0, E v ), 
p=(M, 0, 0, 0), 
k! = (E T . p T s'm6 T , 0, p T cos# T ), 



(2) 



and the following Lorentz invariant variables are 
defined 
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Each subprocess is distinguished by the hadronic 
invariant mass W: W = M for QE, M + < 
W < Went f° r RES. W C ut is an artificial bound- 
ary, and we regard that DIS process occurs in the 
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regions ofW > W cut . We take W cut = 1.4 GeV 
in this report. 

The differential cross section and the spin po- 
larization vector of produced r~ are obtained in 
the laboratory frame as 

da T 
dE T d cos 9 T 
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+ W 2 (p T + E T cos9 T ) 
({E v + E T )p T - 



-^W 5 cos9 T 
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(E V E T +p 2 T )cos6 T )}, 
(5c) 

where Gp is the Fermi constant and k = 



M 



M 2 W /{Q 2 + M 2 W ). 



s= (s x , s y , s z ) is defined in the 



r rest frame in which the z-axis is taken along 
its momentum direction and the y-axis is along 
k x k', the normal of the scattering plane, in the 
laboratory frame. It is normalized as \s\ — 1/2 
for pure spin eigenstates. Wi—i^...^ are structure 
functions defined with the generic decomposition 
of the hadronic tensor, 



Wy, v (p, q) 



' M 2 



W 2 
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p a q 



W 3 



W 5 , (6) 



M 2 * ' 2M 2 
where the totally anti-symmetric tensor e M „ Q/ 3 is 
defined as e i23 = 1- These functions can be es- 
timated for each process, QE, RES, and DIS. See 
explicit forms in Ref. [B]. 
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Figure 1. Contour plot of the DIS cross section 
for the v T N — » t~X process at E v = 10 GeV in 
the laboratory frame. The kinematical boundary 
is shown by a thick curve. The t polarizations 
are shown by arrows. The length of the arrows 
gives the degree of polarization, and the direction 
of the arrows gives that of the r spin in the r rest 
frame. The size of the 100% polarization (P = 1) 
arrow is shown as a reference. 



(i) at a fixed neutrino energy 

In Fig. [TJ the differential cross sections 
da T /dp z dpT, obtained from Eq. ([5]), at the in- 
cident neutrino energy E v = 10 GeV are shown 
as a contour map0, where p z — p T cos 9 T and 
Pt — p T sin 9 T . Only the contours of the DIS cross 
section are plotted to avoid too much complexity. 
Each contour gives the value of the differential 
cross section in the unit of fb/GeV 2 , e.g., the in- 
nermost line is for 4 fb/GeV 2 . The QE process 
contributes along the boundary, and the RES con- 
tributes just inside of the boundary. The contour 
shows that the contributions in the forward scat- 
tering angles in the larger p T side are important. 
In that region, the cross sections of QE and RES 
are also large and comparable to that of DIS. 

The polarization vector s, Eq. ([5|), of t~ is also 
shown in Fig. [TJ The length of each arrow gives 
the degree of polarization (0<P = 2|s|<l)at 
each kinematical point and its orientation gives 
the spin direction in the r rest frame. The pro- 



1 It must be noted that this contour map differs from Fig. 3 
of Ref. where da T /dE T dcos9 T is plotted. 
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Figure 2. Initial flux (left) and v T flux followed 
by — > v T oscillation at Gran Sasso (right). 
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Figure 3. Contour plot of the number of t~ pro- 
duction events in the CNGS experiments. The 
r~ polarization are shown by the arrows. 



duced r~'s have high degree of polarization, but 
their spin directions significantly deviate from the 
massless limit predictions, where all r~'s should 
be purely left-handed. Since s x of Eq. ([5]) turns 
out to be always negative, the spin vector points 
to the direction of the initial neutrino momentum 
axis. Qualitative feature of the results can be un- 
derstood by considering the hclicity amplitudes 
in the center-of-mass (CM) frame of the scat- 
tering particles and the effects of Lorentz boost 
from the CM frame to the laboratory frame. See 
more details in Ref. [5]. Let us stress that these 
features of the polarization of play an impor- 
tant role in the following analysis. 

(ii) in the CNGS experiments 

In the CNGS experiments, beam is produced 
at CERN-PS, which is expected to deliver 4.5 x 
10 19 protons on target (p.o.t.) per year. The 
beam is optimized for v T appearance with (E v ) = 
17 GeV. Figure [2] shows the expected v T flux, 

4>u r (E u )= J2 4>»AE»)xPu^ Vt (E v ), (7) 

£=e, fj, 

at Gran Sasso with the baseline length of L = 
732 km from CERN. Here <j) Ve {E v ) are the ini- 
tial vi fluxes {I = e, n) pQ, see Fig. [2] (left), and 
Pv e — >v r (Ei/) are the vi v T oscillation probabil- 
ities in the three-neutrino model. The fraction 
Vejv^ in the initial fluxes is less than 1%. In our 
analysis, we take the following values for the neu- 



trino oscillation parameters: 

^ m i2.i3 = (8-2 x 1CT 5 , 2.5 x 1CT 3 ) eV 2 , 

sin 2 2012,23,13 = (0.8, 1, 0), 6 M ns = 0°, (8) 

with the constant matter density of p = 3 g/cm 3 . 
Here we assume the so-called normal hierarchy. 

Taking into account the CNGS neutrino flux 
shown in Fig. [2l we show the distributions of 
events and polarization vectors of r~ on the 
p T cos 9 T -p T sin 9 T plane in Fig. [3] The initial neu- 
trino energy is integrated out with the incoming 
v T flux, (t>v T {E v ) of Eq. ([7]), whereas it is fixed 
at 10 GeV in Fig. [TJ The number of r~ pro- 
duction events for all the QE, RES and DIS pro- 
cesses is included in the contour map, where we 
assume 5 years with 4.5 x 10 19 p.o.t. /year of the 
primary proton beam and the 1.65 kton size de- 
tector, which are the current plan of the OPERA 
experiment [3]. Each contour gives a number of 
events per GeV 2 , e.g., the innermost line is for 
7 events/GeV 2 . The contour map shows that 
there are many events around E T = 10 GeV, and 
around 9 T = 5°. As for the polarization vectors, 
the dependence on the energy and the scattering 
angle of t~ is rather smooth as compared to that 
in Fig.[T]because of the integration of the incident 
neutrino energy. However, the direction of the r 
polarization is still non-trivial in the region which 
has many events. 
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3. TAU DECAY DISTRIBUTIONS 

We present our results of the decay particle dis- 
tributions from r~ leptons produced by the CC 
interactions for the CNGS experiments. Main 
feature of our analysis is to deal with the proper 
spin polarization of t~ which is calculated for 
each production phase space, shown in Fig. [3] 
In order to show the effects of the r polariza- 
tion on the decay distributions, we compare the 
results with unpolarized r decays and also with 
completely left-handed r decays. 

The events of the decay distributions for i = 
7r, i are given by 



u 4 



dNi 



,1 / dE v <t>„ T {E v ) 



E thr 



£4 



dcos6 T / dE T - 

J E _ dE T d cos v T 

1 dT; 



Y T dEidfli 



(E T ,9 T ,s{E v ,E T ,9 T )), (9) 



where A is the number of active targets, E™ ax 
is the maximum value of neutrino energy in the 
flux, El hr = m T + m 2 /2M is the threshold energy 
to produce r lepton off a nucleon, and the other 
integral ranges are given by 



E±{E v ,e r ) = (b±Vb 2 

where 




{E v + M) 2 - El cos 2 6 T , 
{E v +M)(ME v +m 2 T /2), 



(12) 



c = m 2 T E 2 v cos 2 9 T + (MEy + m 2 /2) 2 . 



Here, dTi/dEidSli is the decay distribution for i = 
7r, t in the laboratory frame. 

Figure Q] shows the energy distributions of 7r~ 
(thick lines) and i~ (= e _ , (medium-thick 
lines) decayed from produced in the neutrino- 
nucleon CC interactions. We assume the same 
configuration of the experimental setup as Fig. [3l 
i.e., 5 years running with 4.5 x 10 19 p.o.t. per 
year of the primary proton beam and 1.65 kton 




Figure 4. Energy distributions of 7r~ (thick lines) 
and £~(= e~, (jT) (medium-thick lines) in the de- 
cay of r~ produced in the neutrino-nucleon CC 
interactions for the CNGS experiments. Solid, 
dashed, and dotted lines shows the energy dis- 
tributions with the predicted r polarization, un- 
polarized, and purely left-handed cases, respec- 
tively. The estimated number of r~ production 
is also shown by a thin solid line with respect to 
the r energy (E T ). 



size detector for the OPERA experiment [3]- For 
each decay mode, solid lines show the distribu- 
tions from the decay of r~ with the predicted r 
polarization. For comparison, dashed and dotted 
lines show those of unpolarized r and purely left- 
handed r, respectively. The estimated number of 
the production with respect to the r lepton 
energy, E T , is also plotted as a thin solid line. 
The results are calculated by using Eq. ^ with 
100% particle detection efficiency for simplicity. 
For the above parameters, 113 events of t~ are 
produced and 13 (20) of those decay into 7r~ (£~) 
mode. 

The 7T _ and l~ distributions have peak in the 
low energy region, and in this region the polar- 
ization dependence becomes large. The polariza- 
tion dependence is opposite between 7r _ and l~ , 
and is more significant in tt~ than in l~ . In the 
peak region, the polarization dependence affect 
the distribution around 30% (15%) for the 7r~ 
(£~) decay mode. Expected statistics is rather 
small in the current design of the CNGS exper- 
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Figure 5. Azimuthal angle distribution of tt 
(thick lines) and l~ (thin lines). 



iments. However, the likelihood probability of 
each event will be affected significantly by the r 
polarization effects. The characteristic feature of 
our prediction is that the produced t~ is almost 
fully polarized and that it has non-zero trans- 
verse component of the spin vector, namely s x 
of Eq. ([5]). The observed patterns of the ir~ and 
t~ energy distributions in the laboratory frame 
then follow from the energy angular distributions 
in the polarized r rest frame. 

Figure[5]shows the azimuthal angle distribution 
of 7r~ (thick lines) and l~ (thin lines). The az- 
imuthal angle ipi =Vi g is given by d£l\ = dcos B\ dtpi 
in Eq. ([5]), and is measured from the scattering 
plane where ipi = 7r/2 is along the p„ x p T direc- 
tion in the laboratory frame, in which the z-axis 
is taken along the direction of the r momentum. 
Solid lines show the distributions from t~ with 
the predicted t polarization, and dashed lines 
show those from unpolarized r. The results of 
purely left-handed t are the same as those for the 
unpolarized r. Since both unpolarized and purely 
left-handed r~ have zero component of perpen- 
dicular polarization, they give flat azimuthal dis- 
tributions. The azimuthal angle distributions can 
be measured by tracking the trajectory of r lep- 
tons by emulsion detectors in the OPERA experi- 
ment, or by reconstructing the hadronic cascades 
from neutrino-nucleon scattering. As is the case 
of energy distribution, 7r~ and t~ decay mode 



show the opposite feature and polarization de- 
pendence is clearer on 7r~ mode than £~ mode. 
At ip = or \ir\, the dependence of the r polariza- 
tion affects the distribution by about 47% (f6%) 
for the 7T _ (£~ ) decay. Even though the number 
of event is limited, it may be possible to obtain a 
hint of such large asymmetries. 

4. SUMMARY 

In this report, we studied the effects of the spin 
polarization of t~ produced in neutrino-nucleon 
scattering on the subsequent decay distributions, 
for t~ into 7r~ or t~ (= e~, ^r) modes. The spin 
polarization of r _ produced via neutrino were 
briefly reviewed. Taking into account the po- 
larization of produced r _ , we showed the energy 
and azimuthal angle distributions of 7r~ and £~ in 
the laboratory frame, for the experimental setup 
of the CNGS long baseline project, OPERA and 
ICARUS experiments. 

We found that the decay particle distributions 
in the laboratory frame are significantly affected 
by the r~ polarization. Rather strong azimuthal 
asymmetry of 7r~ and £~ about the r momentum 
axis is predicted, which may have observable con- 
sequences even at small statistics experiments. 

Acknowledgements 

K.M. would like to thank M. Sakuda for sup- 
port and encouragement. 

REFERENCES 

1. CNGS project home page, 
http://proj-cngs.web.cern.ch/ 

2. ICARUS Collaboration home page, 
|http://pc nometh4. cern.ch/| 

3. OPERA Collaboration home page, 
http://operaweb.web.cern.ch/ 

4. K. Hagiwara, K. Mawatari and H. Yokoya, 
Nucl. Phys. Proc. Suppl. 139 (2005) 140. 

5. K. Hagiwara, K. Mawatari and H. Yokoya, 
Nucl. Phys. B 668 (2003) 364; Erratum-ibid. 
B 701 (2004) 405. 

6. M. Aoki, K. Hagiwara, K. Mawatari and 
H. Yokoya, Nucl. Phys. B 727 (2005) 163. 



